Determination of the magnetic domain size in the ferromagnetic superconductor 
UGe 2 by three dimensional neutron depolarization 
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Three dimensional neutron depolarization measurements have been carried out on single- 
crystalline UGe2 between 4 K and 80 K in order to determine the average ferromagnetic domain 
size d. It is found that below Tc = 52 K uniaxial ferromagnetic domains are formed with an 
estimated magnetic domain size of d ~ 4 — 5 /im. 

PACS numbers: 61.12.-q, 75.50.Cc, 75.60.Ch 
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I. INTRODUCTION 



Recently, the compound UGe2 has attracted much at- 
tention because superconductivity was found to coexist 
with ferromagnctism^^ Until this discovery, only su- 
perconducting compounds exhibiting antiferromagnetic 
order had been known like DyMogSg, GdMoeSg, and 
TbMogSgi^^ Coexistence of antifcrromagnetism and su- 
perconductivity was also found in the heavy fermion com- 
pounds like Celn 3 , CePd 2 Si 2 , and UPd 2 Al 3 &La In these 
cases superconductivity and antifcrromagnetism appear 
simultaneously, because the Cooper pairs are insensitive 
to the internal fields arising from the antiferromagnetic 
ordering when the superconducting coherence length £ is 
much larger than the periodicity of the static antiferro- 
magnetic ordered structure. However, in a ferromagnetic 
structure we expect that the internal fields do not can- 
cel out on the length scale of £ and therefore have their 
influence on the Cooper pairs. That ferromagnetic or- 
der excludes superconductivity, is nicely demonstrated 
in ErRh 4 B42iA& and HoMo 6 Sa±I, where standard BCS 
singlet-type superconductivity is suppressed when ferro- 
magnetic order sets in. Otherwise, if one would consider 
unconventional spin-triplet superconductivity mediated 
by ferromagnetic spin fluctuations, the pairing is rela- 
tively insensitive to a local magnetic field and can there- 
fore coexist with ferromagnetic order. On the other hand, 
when the ferromagnetic domain size d is much smaller 
than the superconducting coherence length £, one effec- 
tively has no internal magnetic field. 

The coherence length £ for UGc 2 is cstimatedi*i£ to be 
130-200 A. Interestingly, Nishioka et ajpiii4, considering 
jumps in the magnetization at regular intervals of mag- 
netic field and at very low temperatures, estimated the 
ferromagnetic domain size d to be of the order of 40 A, by 
attributing the jumps to field-tuned resonant tunnelling 
between quantum spin states. Since d would be several 
times smaller than £, it was proposed that the ferromag- 
netism can be cancelled out on the scale of the coherence 
length of the Cooper pairs. This would imply that the 



pairing mechanism for superconductivity might be of the 
singlet-type after all. 

UGe 2 crystallizes in the orthorhombic ZrGa 2 crystal 
structure (space group Cmmm)i^ with unit cell dimen- 
sions a = 4.036 A, b = 14.928 A, and c = 4.116 A, con- 
taining 4 formula units. Ferromagnetic order sets in at 
Tc = 52 K. The saturated magnetic moment at ambient 
pressure is 1.4 /xb/U, directed along the a axis. Magnetic 
measurements indicate a very strong magnetocrystallinc 
anisotropy^ Superconductivity is found only in a lim- 
ited pressure range between 10 and 16 kbar with a max- 
imum transition temperature of T c ss 0.7 K. In this pres- 
sure range, the magnetic moment is still 1 //b/U. Within 
the ferromagnetic phase, a second transition occurs at 
Tx ~ 25 K at ambient pressure, below which the mag- 
netic moment is enhanced. Therefore the temperature 
region from Tc to Tx was named the weakly polarized 
phase, while the lower temperature region T < T\ was 
called the strongly polarized phased 

In this paper we report on three dimensional neu- 
tron depolarization measurements performed on single- 
crystalline UGe 2 at ambient pressure between 4 K and 
80 K. Our principal aim was to determine the ferromag- 
netic domain size d in UGe 2 and compare the value to 
the size of ~ 40 A estimated by Nishioka et alm& on the 
basis of the hypothesis of field-tuned resonant tunnelling 
between spin quantum states. Since the neutron is a very 
sensitive probe to local magnetic fields, neutron depolar- 
ization is an excellent technique to measure the average 
domain size and the domain wall width. 



II. EXPERIMENT 
A. Experimental 

The measurements were performed on the Poly Axis 
Neutron Depolarization Analyzer (PANDA) at the In- 
terfaculty Reactor Institute (IRI) of the Delft University 
of Technology. The used neutron wave length was 2.03 A 
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which corresponds to a velocity of 1949 m/s. 

The neutron depolarization measurements on UGe2 
were performed on a single-crystalline sample with di- 
mensions a x b x c = 4.0 x 0.440 x 3.0 mm 3 . The b 
axis was oriented along the transmitted neutron beam 
(x) with a transmission length L and the easy axis for 
magnetization a along the vertical axis (z) within the 
plate of the sample. The crystal has been grown from 
a polycrystalline ingot using a Czochralski tri-arc tech- 
nique. No subsequent heat treatment was given to the 
crystal. The illuminated area was a rectangle with di- 
mensions y x z = 1 x 2 mm 2 centered at the middle of 
the sample. 

The measurements in zero field were performed during 
a temperature sweep from 2 K up to 80 K and down to 
2 K with a low sweep rate of 10 K/hr. The measurements 
in non-zero field (4 and 8 mT) were done during a similar 
temperature sweep with a sweep rate of 25 K/hr. The 
sample was first zero- field cooled, whereafter the field 
was switched on at the start of the measurements. The 
subsequent measurements were performed during heating 
and cooling in a constant field. 



B. Neutron depolarization 

The neutron depolarization (ND) technique is based on 
the loss of polarization of a polarized neutron beam af- 
ter transmission through a (ferro)magnetic sample. Each 
neutron undergoes only a series of consecutive rotations 
on its passage through the (ferro)magnetic domains in 
the sample. It is important to note that the beam cross 
section covers a huge number of domains, which results 
in an averaging over the magnetic structure of the whole 
illuminated sample volume. This averaging causes a loss 
of polarization, which is determined by the mean domain 
size and the mean direction cosines of the domains. The 
rotation of the polarization during transmission probes 
the average magnetization. 

The 3x3 depolarization matrix D in a ND experiment 
expresses the relation between the polarization vector be- 
fore (P°) and after (P 1 ) transmission through the sample 
(P 1 = DP ). The polarization of the neutrons is cre- 
ated and analyzed by magnetic multilayer polarization 
mirrors. In order to obtain the complete matrix D, one 
polarization rotator is placed before the sample and an- 
other one right after the sample. Each rotator provides 
the possibility to turn the polarization vector parallel or 
antiparallel to the coordinate axes x, y, and z. The re- 
sultant neutron intensity is finally detected by a 3 He de- 
tector. The polarization rotators enable us to measure 
any matrix element Dij with the aid of the intensity of 
the unpolarizcd beam 1$: 



The matrix element £),-,- is then calculated according to 



h 



(1) 



where ly is the intensity for P° along i and P 1 along j. 



Da = 



1 Is-h 



P> Is 



(2) 



where Po is the degree of polarization in the absence of 
a sample. In our case we have Po = 0.965. 
We now introduce the correlation matrix ay : 

^dx'AB i (x,y,z)AB j (x',y,z)), (3) 

where AB(r) = B{f) — (B) is the variation of the mag- 
netic induction, (. . .) denotes the spatial average over the 
sample volume and where the integral is over the neutron 
transmission length L through the sample. Assuming 
= for i ^ j we define the correlation function £ as 



(4) 



With these two quantities it can be shown that if there 
is no macroscopic magnetization ((B) = 0) the depolar- 
ization matrix is diagonal and under the assumption of 
atj = for i ^ j given by 18 i 19 i 20 



Di 



x,y,z. 



(5) 



where 7 = 1.83 x 10 s T is the gyromagnetic ratio 
of the neutron and v its velocity. 

Intrinsic anisotropy is the depolarization phenomenon 
that for magnetically isotropic media the depolarization 
depends on the orientation of the polarization vector with 
respect to the propagation direction of the neutron beam. 
The origin of this intrinsic anisotropy is the demagneti- 
zation fields around magnetized volumes in the sample. 
In the following we will assume that the demagnetization 
fields are negligible for needle shaped magnetic domains. 

We now discuss the case (B) 7^ 0. When the sample 
shows a net magnetization, the polarization vector will 
rotate in a plane perpendicular to the magnetization di- 
rection. If the sample shape gives rise to stray fields, the 
rotation angle <f> is related to the net magnetization (M) 
by 



= r]-Ln (M) = T]^Ln M s (m) (6) 
v v 

where 77 is a geometrically factor given in Eq. EH for 
a rectangular shaped sample and (m) = M/M$ the re- 
duced sample magnetization in terms of the saturation 
magnetization Ms = M$(T). If the mean magnetic in- 
duction (B) in the sample is oriented along the z-axis, the 
depolarization matrix is, for <j) 3> (7/w) \ot X x — a yy \ L/2 
(the weak damping limit), given bj i i2 4 i2i2£ 
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FIG. 1: The diagonal elements of the measured depolariza- 
tion matrix D for increasing and decreasing temperature for 
UGe2- All other elements of the depolarization matrix are 
zero within the experimental uncertainty. For D xx and D yy 
the experimental uncertainty is within the symbol size. 



With the net magnetization along the z-axis, the ro- 
tation angle <f> of the beam polarization is obtained from 
the measurements by 



arctan 



and £ is calculated with 



D, 



D, 



D 



HI) 



i = -v 2 ln{detL>}/2 7 2 £. 



(8) 



(9) 



As mentioned earlier, ND also provides information 
about the mean-square direction cosines of the magnetic 
induction vector in the (ferro)magnetic domains. These 
are directly given by the quantities 7; = a«/£, where 
i = x,y, z, and can be estimated from the measurements 
by 



7i = 1 - 2 In Da/ In {det D} 



(10) 



This equation is only valid for those directions which 
show no net rotation of the beam polarization. 



III. RESULTS 

A. Measurements in zero field 

In Fig. 2]wc show the diagonal elements of the depo- 
larization matrix for UGe2, measured in zero magnetic 
field. All off-diagonal elements are zero within the ex- 
perimental uncertainty in the studied temperature range. 
The measurements for increasing temperature are quali- 
tatively the same as those for decreasing temperature, as 
expected. 

The Curie temperature of Tq = 52 K is clearly indi- 
cated in Fig. nby the kink in D xx and D yy . Note that 
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FIG. 2: The determinant of the measured depolarization ma- 
trix det D and the rotation angle (f> of the beam polarization 
after passage through the sample of UGe2 in 4 mT and 8 mT 
for increasing and decreasing temperature. 



Dyy below Tc indicates that there is no intrinsic 



anisotropy and hence that the magnetic domains produce 
virtually no stray fields. Furthermore, D zz w 1 indicates 
that all moments arc oriented along the a axis. 



B. Measurements in small field 

In Fig. |3 we show the determinant of the depolariza- 
tion matrix det D and the rotation angle cj> after passage 
through the sample for measurements in fields of respec- 
tively 4 and 8 mT (after zero- field cooling). The data 
of (f> have been corrected by subtracting the mean value 
above Tc, since this rotation is merely due to the applied 
field. 

At low temperatures the magnetic fields (applied after 
zero field cooling) are too small to fully align the mag- 
netic domains. Therefore, the measurements for increas- 
ing and decreasing temperature do not yield the same 
results. Whereas for increasing temperature the rota- 
tion shows an increase, for decreasing temperature the 
data represent a monotonous magnetization curve, as ex- 
pected for a field-cooled ferromagnct. Close to Tc there 
is no difference between field cooling or field warming. 

Fig. [21 shows that for 4 mT the depolarization is at 
the same level as for mT. Above T\, however, ex- 
tra depolarization occurs. This means the system gets 
more inhomogeneous, i.e. the domains grow and the 
magnetic correlation length increases (£ in Eq. 0), lead- 
ing to extra depolarization. Close to Tc the depolariza- 
tion disappears, because the magnetic moment decreases 
sharply. For decreasing temperature the determinant has 
the same shape as in the case of mT. At 8 mT the de- 
terminant is already reduced below Tx , indicating larger 
domains. 

Again, the Curie temperature of Tc = 52 K is clearly 
indicated by the kink in det-D and <j). Also note the 
abrupt increase in <j) around Tx ~ 20 K. Evidently 
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the system passes, with increasing temperature, from a 
strongly polarized phase to a weakly polarized phase, as 
reported earlierpii 



IV. DISCUSSION 



Model and results 



The measurements confirm that UGC2 is a highly 
anisotropic uniaxial ferromagnet. Further, the magnetic 
domains are long compared to their (average) width, be- 
cause D xx = Dy V indicates relatively weak stray fields 
produced by the magnetic domains. This allows us to 
assume B{r) = fj,oM(r) inside the domains. In order to 
analyze the data we consider a model where the sample 
is split into -/V long needles along the a axis with a fixed 
width A and a magnetic induction B$ = fJ-oMs along 
the a axis. With (iV|J the number of domains with 
a magnetic induction pointing upwards (downwards), we 
can define the reduced macroscopic magnetization of the 
sample, pointing along the z-direction, as 



iV T - N± _ (B z ) 
iV T + Ni B s 



(11) 



Each needle will have magnetic induction f or [ with 
probability pf = (1 + (m z ))/2 and p^ = (1 — (m z })/2, 
respectively. The polarized neutron beam traversing the 
sample will therefore see a binomial distribution of f and 
|, which results in a depolarization matrix D with ele- 
ments 



D, 



^A(l-(ra«) 2 ) ( jB S L 

\ v ' ' cos I \m z ) I , 



D- 



2« 2 V v *' ) sin \ m z) 

\ v J 



(12) 



and all other elements equal to 0. (Note that, since we 
have not taken into account the macroscopic stray fields, 
the angle <j) should be corrected by the factor of r) (Eq. |fjj 
before calculating (m z ) in Eg. 1121) 

Within this binomial distribution model it is easy to 
show that for the case (m z ) = the average ferromag- 
netic domain size d is equal to 2 A. Given a domain 
wall (i.e. two adjacent needles with an opposite mag- 
netic induction), the probability of forming a domain 
of n needles is (^) and the average is calculated by 



= 2. When a field is applied, we have to 



distinguish between a domain (with size d-\ ) in which the 
magnetic induction is parallel to the field and a domain 
(with size d^) with opposite induction. The probability 
of forming a domain of size n is p r ?~ 1 pi = p*~ (1 — pj) 
which leads to d T /A = 1/(1 -p T ) = 2/(1 - (m z )). Simi- 
larly, d ; /A = 2/(1 + (m,». 

In order to estimate the domain wall thickness 8 we 
assume that m z changes sinusoidally from +1 to -1 over 
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FIG. 3: Calculated values of the reduced macroscopic mag- 
netization (m-) and the average ferromagnetic domain sizes 
with magnetic induction parallel (dj) or antiparallel (djj to 
the applied magnetic field of 4 mT and 8 mT in UGe2 for 
increasing and decreasing temperature. 



a distance S in the form of a Bloch wall. The consequence 
is that D zz is slightly less than 1 in the ordered state. For 
such a domain wall it is straightforward to show that the 
domain wall thickness 5 can be estimated by 



(13) 



which can be measured directly by Eq. 1101 

For the values of Bs needed in Eq. ^1 we use the ex- 
perimental magnetic moment of Ref^i which we convert 
to magnetic induction, remembering there are 4 formula 
units per unit cell. For the value of r\ in Eq. I^Jwe take 
i] = 0.6. 

From Fig. [I] it is clear that the data for increasing and 
decreasing temperature give slightly different results for 
the ferromagnetic domain size d in zero magnetic field. 
The values found for d = 2 A are 5.f (2) /jm when cooling 
down slowly and 4.4(1) /xm when heating up after fast 
cooling. Both values are independent of temperature. 
These values indicate the domain size perpendicular to 
the a axis (along the b axis). Along the a axis which we 
assume the domain size is much larger. 

The magnetic domain wall thickness S, divided by the 
magnetic domain size d, is calculated with Eq. 1131 from 
the experimental data in Fig. ^ and amounts to S/d = 
0.047(23), independent of temperature. This gives 5 = 
0.22 /im. The size of the domain wall thickness is thus 
found to be only a minor fraction of the domain size. 

Analysis of the data in a small magnetic field (Fig. 0) 
with Eq. 1121 gives the results shown in Fig.[3]and tabled 
For 4 mT and increasing temperature (after zero field 
cooling), the reduced magnetization (m z ) remains equal 
to up to Tx ~ 20 K. As a consequence d± is equal to 
d-f and of the same order of the zero field values. Above 
Tx however the system gets magnetically soft and (m z ) 
starts to increase linearly towards ~ 0.7. Domain walls 
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TABLE I: Ferromagnetic domain sizes in UGe2 for increas- 
ing temperature after zero field cooling (ZFC) and decreasing 
temperature in field (FC). The sizes of the domains with mag- 
netization parallel to the applied magnetic field is denoted 
by d| and the domains with antiparallel magnetization by 
dj. Below Tx the domain sizes are temperature independent. 
Above Tx the domains grow. The values shown are at a few 
Kelvin below Tc- 
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are expelled above Tx, since increases much faster 
than d\. (Note the logarithmic vertical scale.) While c?| 
gets of the order of 100 /im, di only reaches 25 /im. When 
the domains grow in width, at a certain moment it is no 
longer allowed to assume B{r) = noM(r), because stray 
fields produced by the domains have to be taken into 
account. The model therefore is no longer appropriate 
close to Tc- 

For field cooling in 4 mT, the system has (m z ) = 
0.668(1) for the whole temperature range below Tc- The 
values of the domain size are shown in table 

When after zero field cooling a field of 8 mT is turned 
on, the sample does get a macroscopic magnetization, 
in contrast to the case of 4 mT. Up to Tx ~ 20 K the 
reduced magnetization (m z ) = 0.221(2) is independent of 
temperature. Then (m z ) starts to increase up to 0.718(3) 
around 30 K and is constant afterwards up to Tc- When 
cooling down in 8 mT, (m z ) = 0.708(1) over the whole 
temperature range below Tc- 

For field warming after zero field cooling, the calcula- 
tion of the domain sizes yields unexpected temperature 
dependencies of domain sizes above Tx ■ As can be seen 
in Fig. |3 according to the model the domain sizes grow 
above Tx to decrease in size again at higher temperature. 
Clearly there is another source of depolarization, not ac- 
counted for by the model. Since the field is strong enough 
to penetrate the sample, additional depolarization arises 
from an inhomogeneous magnetic domain structure. 



B. Discussion 

If the domain width becomes relatively large compared 
to its length, stray fields become important and the sim- 
ple model assuming B(r) = [ioM(r) is no longer valid. 
Calculation of the mean-square direction cosine along the 
z direction, 7,, with Eg. 1101 indeed shows a decrease from 
unity above Tx, indicating that the magnetic induction 



B is not along the a axis throughout the sample. The 
model can of course be improved if we no longer assume 
a length/ width ratio of infinity (no stray fields) for the 
domains. The simple model together with our measure- 
ments, however, do show that the magnetic domain sizes 
in zero field are a few /im and that by applying small 
fields the domains grow. Our measurements therefore 
indicate that the domain sizes in UGe2 at ambient pres- 
sure and down to 2 K are certainly larger than the 40 A 
predicted by Nishioka et alM^ 

In Fig. ^ it is shown that D zz is less than unity be- 
low Tc- This can be caused by the domain walls, but 
can also be accounted for by a misalignment. A simple 
calculation shows that a misalignment of 8° would fully 
account for the values of D zz below Tc- The stated value 
of S = 0.22 /xm (or S/d = 0.047(23)) should therefore be 
regarded as an upper limit. 

From the above considerations we conclude that the 
domain structure of UGe2 behaves like in a conventional 
ferromagnet. The magnetic domain size largely exceeds 
the superconducting correlation length of the Cooper 
pair. The magnetic domain boundaries can therefore 
only give secondary effects on the superconducting or- 
der. 



V. CONCLUSION 

The ferromagnetic domain sizes of UGe2 was studied 
by means of three dimensional neutron depolarization 
at ambient pressure. We conclude that the existence 
of field-tuned resonant tunneling between spin quantum 
statesi^i^ is highly unlikely. The requirement of this 
model is a ferromagnetic domain size of 40 A while our 
measurements indicate a size a factor of 1000 larger. 
The observed jumps in the magnetization should be at- 
tributed to a Barkhauscn effect as discussed by Lhotcl 
et alS^k. The superconductivity therefore exists within a 
single ferromagnetic domain. The domain walls are not 
expected to strongly affect the bulk Cooper pair wave 
function, as suggested by Nishioka et alm&l^, since the 
domain wall is less than a few percent of the average 
domain size. 



APPENDIX A: EFFECT OF STRAY FIELDS 
INDUCED BY A HOMOGENEOUSLY 
MAGNETIZED SAMPLE 

In this appendix we will calculate the magnetic induc- 
tion B generated by a uniformly magnetized sample with 
length I, width w, and thickness t (Fig.0J. Moreover, an- 
alytical expressions will be given for the line integrals of 
B along the path of a neutron. The center of the sample 
is taken as the origin of the reference frame. 
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FIG. 4: Schematic layout of a homogeneously magnetized 
rectangular sample. 



Our starting point is the Biot-Savart law: 

Mo f VxMxf 



B(x,y,z) = — I 3 dS 



-in 



dr 



(Al) 

where /io = 47r x 10 -7 H/m, M is the magnetization, 
n the unit vector perpendicular to the surface S of the 
sample, r the vector pointing from the surface S to the 
point (x, y, z), and r the volume of the sample. Since 
the sample has a homogeneous magnetization, the second 
term vanishes and, with rfi = M/M, we have 



g _ Mo A/ f mxnx f 

47T 



(A2) 



A straightforward but tedious calculation shows that 
B x = - ^ eie 2 £3ln(<S- {2y + e 2 w)) 

ei,e2,e3 — ±1 

B y = - J! ei£2e 3 ln(5- {2x + at)) (A3) 

ei,e2,e3=±l 

2z + e 3 l 



D 



where 



eie 2 £3 arctan 



£l,£2,e3=±l 



(2a + eit) (2y + e 2 u;) 
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S = yj (2a + eii) 2 + (2y + e 2 w) 2 + S + esO 2 - (A4) 

The rotation of the polarization of a neutron beam 
depends on the line integral of the magnetic field along 
the neutron path. From the Larmor equation 4zP(t) = 
jP(t)xB(t), or equivalents ^P(x) = (7/w) P{x)xB(x) 
where 7 = 1.83 x 10 s s _1 T _1 the gyromagnetic ratio of 
the neutron and v its velocity, we get the general solution 



P(x,y,z) 



I - / B(x',y,z)dx' ) } P(-oo,y,z) 



(A5) 



where we have defined the magnetic field tensor B as 
f B z —B v 



B(x,y,z) 



\ 



—B z 

By —B X 



B x 




(A6) 



Thus, in order to calculate the rotation of the neutron 
beam polarization due to a homogeneously magnetized 
sample, the following line integrals arc required: 



X(y,, 



B x (x',y,z)dx' = 



Y (y, z ) = / B y (x',y,z)dx' 

J — OO 

tfioM 



-in 



^2 e 2 e 3 In ({2y + e 2 w) 2 + {2z + e 3 lf 



z (y, z ) = I B z (x' ,y 7 z)dx' 

00 

^2 6263 arc ta- n 



t(i M x - / 2z + e 3 Z 



2tt 



2y + e 2 w 



(A7) 



For completeness we also give the line integrals in the 
case the neutron beam is along the z-direction: 



X'(x,y)= / B x (x,y,z')dz' = 

J — OO 

/OO 
B y (x,y,z')dz' = 
-OO 
/* OO 

Z'(x,y)= B z (x, y, z')dz' = fi M I (A8) 



From Eq. IA5I and the above line integrals we get for 
the final polarization P(oo, y, z) = D(y, z) x P(— 00, y, z) 
with D(y, z) equal to 

(Ecosav 7 ^ ZVZsinaVT, -r v / Esinav / E \ 

-ZVHsinay/Z Z 2 cos aVS + Y 2 YZ (l — cos a Ve) 
Fv^sinav^ yZ(l-cosa V / E) F 2 cos av 7 ^ + Z 2 / 

(A9) 

where z) = Y 2 (y, z) + Z 2 (y, z) and a = 7/f . 

Now Eq. IA9I relates the initial polarization to the fi- 
nal polarization for a beam passing through the sample 
at (y,z). For a neutron beam with finite cross section 

the matrix D should be integrated over the beam cross 
section. If the integration is symmetric relative to the 
origin, then we can make use of the fact that 



B y {eix, e 2 y, £3z) = e 2 e 3 B y {x, y, z), 



(A10) 



where ei,e 2 ,£3 = ±1. This means that B y (x,y, z) and 
therefore Y(y, z) are antisymmetric with respect to y and 
z. Therefore, from Ecr. lA9l we only have to integrate the 
matrix 



E cos a\/E ZVE sin a\/E 



£ -Zv / Esinav / E Z 2 cos aV^ + Y 2 

Y 2 cos a\/E + Z 2 J 

(All) 

over the cross section of the neutron beam. 

An infinitely narrow neutron beam passing exactly 
through the middle of the sample will only have its 
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polarization rotated by B z (x,0,0) since B x (x,0,0) and 
B y (x, 0, 0) vanish. As long as Y 2 (y, z) is small compared 
to Z 2 (y, z), which is valid if (y, z) is sufficiently far from 
the edges, Eg. lAlil is a pure rotation matrix. 

It is now possible to calculate the magnetization of 
the sample from the measured rotation angle (Eq. |HJ. If 
no stray fields were present, the rotation angle would be 
given by tfigMj /v. However, the stray fields reduce the 
rotation angle to r y/vZ(y, z) with Z(y, z) given in Ea. lA7l 
We can therefore define the geometrical factor rj as 



4> = rjfj.oMt'y/v 



(A12) 



Since 77(0, 0) is a saddle point (77(0, z) has a local maxi- 
mum and 77(2/, 0) a local minimum) an average over the 
cross section of the neutron beam, centered around the 
middle of the sample, will yield a result very close to the 
value of 7/(0,0), which is given by 



2 / I 

— arctan — 

7T \W 



(A14) 



where 77(2/, z) is given by 



V{y, z ) = 7T~ / . £ 2£ 3 arctan — 

£2, £3 



(A13) 



* Electronic address: |l5 . SakaryaQiri .tudelft ,nl| 
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